arXiv:1505.08069v2 [cs.IT] 5Jun2015 


1 


Robust Design of Transmit Waveform and Receive 
Filter For Colocated MIMO Radar 

Wei Zhu and Jun Tang 


Abstract —We consider the problem of angle-robust joint trans¬ 
mit waveform and receive filter design for colocated Multiple- 
Input Multiple-Output (MIMO) radar, in the presence of signal- 
dependent interferences. The design problem is cast as a max- 
min optimization problem to maximize the worst-case output 
signal-to-interference-plus-noise-ratio (SINR) with respect to the 
unknown angle of the target of interest. Based on rank-one 
relaxation and semi-definite programming (SDP) representation 
of a nonnegative trigonometric polynomial, a cyclic optimization 
algorithm is proposed to tackle this problem. The effectiveness 
of the proposed method is illustrated via numerical examples. 

Index Terms —MIMO radar, optimization, receive filter, robust 
design, waveform design. 

I. Introduction 

Due to many advantages over conventional phased-array 
radar [1]—[4], multiple-input-multiple-output (MIMO) radar 
has been widely studied over the last decade. For both colo¬ 
cated [2] and distributed MIMO radar [5], one of the most 
crucial problems is how to design probing signals properly. 
Existing design approaches can mainly be classified into five 
categories according to the criteria adopted: 1) optimizing 
the radar ambiguity function [6], [7]; 2) matching a desired 
beam-pattern [8]—[11]; 3) optimizing the detection or estima¬ 
tion performance based on information theory [12]—[15]; 4) 
optimizing an estimation-oriented lower bound (e.g., Cramer- 
Rao bound [16] and Reuven-Messer bound [17]) and 5) joint 
transmit waveform and receive filter design to maximize the 
signal-to-interference-plus-noise-ratio (SINR) [18]—[21]. 

This letter focuses on the last design approach for colocated 
MIMO radar. In this design framework, joint transmit and 
receive beamforming is investigated in [19] for an active array 
in the presence of signal-dependent interference. A sequential 
optimization algorithm is proposed to maximize the output 
SINR. In [20], joint transmit waveform and receive filter 
design is considered under the constant modulus and similarity 
constraint. Both works rely on exact knowledge of target 
and interferences. Indeed, the angle and INR of interferences 
can be obtained from knowledge-aided methods or estimated 
through previous scans of the space in high INR cases [22], 
[23], The known target angle assumption can be applied to 
the confirmation of an initial detection at some angle bin [24]. 
However, there are other situations where the target angle is 
unknown (e.g., weak target embedded in strong interferences), 
and the SINR should be averagely optimized over the uncertain 
area to avoid beampattern loss [4]. Hence, angular-robust 
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design must be considered and the robust design can also 
be used as an initial step for cognitive detection. In previous 
works [24], [25], robust waveform design has been studied for 
interpulse (or intrapulse) coding in radar by considering the 
unknown Doppler shift of target. Motivated by these works, 
we consider the problem of angular-robust design for colocated 
MIMO radar in the presence of signal-dependent interferences, 
which are induced by the interaction of transmit waveform 
with unwanted scatters. Based on the SINR criterion, transmit 
waveform and receive filter are jointly optimized to maximize 
the worst-case output SINR. Since the resulting problem is 
non-convex, cyclic optimization [26] and semi-definite relax¬ 
ation (SDR) [27] are used to solve it. Although the cyclic 
optimization converges to a locally optimal solution, it still 
can yield a good enough solution with higher worst-case SINR 
than the non-robust design, as illustrated in section IV. This 
is fundamentally different from the optimization problem in 
parameter estimation, in which the local convergence may 
significantly deteriorate the accuracy of estimation. SDR is a 
powerful approximation technique to solve a host of difficult 
non-convex problems with rank constraints. It is commonly 
used in radar signal processing problems, e.g., [22]-[25], 

Notations: Matrices are denoted by bold capital letters, and 
vectors by bold lowercase letters. (-) T , (-) c and {-) H denote 
the transpose, conjugate and conjugate transpose, respectively. 
|| • || denotes Euclidean norm. 0 denotes Rronecker product. 
I[ means Lx L identity matrix. R and C denotes the sets of 
all real numbers and complex numbers, respectively. <5(-) rep¬ 
resents Rronecker delta function, vec (•) denotes vectorization 
operator. Re{-} denotes the real part of the argument. 

II. Problem Formulation 

Consider a colocated MIMO radar system equipped with 
Nt transmitters and Nr receivers. Both the transmit and re¬ 
ceive arrays are assumed to be uniform linear arrays with half- 
a-wavelength element-separation. Let S £ C Nt x N denote 
the transmitted waveform matrix, where N is the number of 
samples in the duration of the transmitted waveform. For a 
particular range cell of interest, the received waveform matrix 

Y £ cVrxiv f rom ]\[ R receivers is corrupted by K signal- 
dependent interferences (e.g., other targets in a multi-target 
scenario [23]) from adjacent range cells with the additional 
noise, and is modeled as 

K 

Y = a 0 a r ( 9 0 )a t (0 o ) T S + ^ a k a r ( 9 k )a t ( 9 k ) T SJ rk +N 

k= 1 

where 
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• (id and a k are the complex amplitudes of the target and 
the fc-th interference source, respectively. 

• 6 *o and O k are the direction of arrivals (DOA) of the target 
and the fc-th interference source, respectively. 

. a r (0) G is the receive steering vector defined by 

a r (9) = [l, e f 7rsin ( 0 ), ■ • • >e MN R - i)sin(@)j' r _ 

. a t (9) G C Ntx1 is the transmit steering vector defined 
by a t (0) = [l,e J7rsin ( 0 ), • • • , e MNr- i)sin( 0 )] T 

• Jr k ,r k G {-N + 1, • • • , -1,0,1, • • • ,N - 1} is an N- 
by-N shift matrix with (Zi.,Z 2 )-th element J r (li,l 2 ) — 
<5(Zi — l 2 — r). r k is the range cell index of the fc-th 
interference source relative to the range cell of interest. 

• N is spatially and temporally white circularly symmetric 
complex Gaussian noise with mean zero and variance cr 2 . 

Let y = vec(Y), s = vec (S) and n = vec (IV). The 
vectorization form of the measurement model is given by 

K 

y = a o A(0 o )s+ '^2a k B(0 k )s + n (1) 

k=t 


where A(0 o ) = I N ® [a r ( 0 o )a t ( d 0 ) T ] and B(0 k ) = Jf k <8> 
[a r ( 0 k )a t [0 k ) T ]. The SINR at the output of the receive filter 

W G C N RNxl j s gj ven by 


X(s,w,0 Q ) 


SNR|ut ff A(6»o)s|“ 

w H Hi(s)w + w H w 


( 2 ) 


where E/(s) = X^^=tINR kB(0 k )ss H B(0 k ) H with the 
signal-to-noise-ratio (SNR) of target and the interference-to- 
noise-ratio (INR) of fc-th interference defined as SNR = 
E{|a 0 | 2 } /c t 2 and INRfc = E{|afc| 2 } /a 2 , respectively. 

We assume that the angle and the INR of the interferences 
are all known or previously estimated, as with prior works 
[19], [20], [23]. We assume that the angle of target under test 
is known to lie in an angular sector 12 = [Oq — A 9, Oq + A 0 ] 
centred around 0(<, where AO indicates the level of angular 
uncertainty. The goal is to maximize the worst-case SINR 
to improve the worst-case detection performance, under the 
waveform energy constraint ||s || 2 = E. Therefore, the robust 
design of transmit waveform and receive filter can be formu¬ 
lated as the following max-min problem: 


max min X(s,w,0q ) subject to ||s || 2 = E (3) 
s,w o 0 en 

Note that for the case of known target angle, (3) reduces to 
the optimization problem in [19]. 


III. Max-Min Robust Design Algorithm 

In this section, we shall present our algorithm to solve 
the problem (3). To begin with, we make some mathematical 
transformations to the objective function of the optimization 
problem. Define W G i^n r xn suc b ^at w = vec(VL). 
Let p(v) = [1, e jV , • • • , e^^ L ~ 1 ' l ] T with L = Nr + Nr — 1 
and v = 7rsin(#o)- Let H = [Hi, Hj, ■ ■ ■ ,H-^ R ] T where 
H G H k G R NtxL , k = 1 ,■■■ ,N r , and the 

(to, n)-th element of H k is defined by H k (m,n) = S(n — 
m — k + 1). Then, one can easily show that a, (0o)®a t (Oq) = 


Hp(v). According to the property of Kronecker products that 
vec (CXD) = (D t ® C)vec(X), we can show that 

w H A(0 o )s = (A(Oo) T vec (W*)) T s (4) 

= vec(a t (0 o )a r (0 o ) T W*) T s (5) 

= {(W H ® I NT )vec(a t ( 0 o )a r (0 o ) T )) T s ( 6 ) 

= s T (W H ®I NT )(a r (0o)®at(0o)) (7) 

= vec(SW H ) T Hp{v) ( 8 ) 

Then, it follows that \w H A(0 q)s\ 2 = p(v) H Gp{v) where 

G = H H vec(SW H )\ec(SW H ) T H (9) 

Let S = [si,S 2 ,--- An] and W = [wi,w 2 ,--- , 

Using SW H = J2n= 1 s nWn' we can write 

N N 

vec (SW H ) = ^ vec (s n w R ) = ^ ut* ® s n (10) 

n= 1 n= 1 

Define X = ss H and V = ww H with X G <C n tNxn t n anc j 
V g C NrNxNrN . Partition X and V into a N-by-N block 
matrix with (m, ri 2 )-th block denoted by X[ ni „ 2 ] G C ,VtX,Vt 
and Vj ni)Tl2 ] G C NrxNr , then it follows that 

G(X,V) = H H ( Vk,« al ®*[ ni .n 2] ) ff (11) 

Xl<rn,n 2 <IV ' 


where we use the notation G(X, V) to emphasize G 
as a function of X and V. Moreover, using w H w = 
tr(V)tr (X)/E and w I1 'i: I (s)w = tr(£,(V)X) with 
S/(V) = j2 k =i mR kB(0k) H VB(0 k ), it is easy to find 
that the denominator of (2) can be re-written as tr((S/(V) + 
U ^'> I Nt n)X) . Consequently, problem (3) can be recast as 


max min 
x,v vex 


subject to 


p{v) H G{X,V)p[v) 
tr ((s ; (y) + fii„)i) 
tr (X) = E, X > 0 , V A 0 
rank (X) = 1, rank (V) = 1 


where X = [uc — A v, vc + Av\ is the corresponding uncertain 
range of v after parameter transformation. 


A. Optimization with respect to X and V 

Since the rank constraint in (12) is non-convex, we adopt 
the commonly-used SDR technique [27] to obtain a relaxed 
problem by dropping the rank-one constraint in ( 12 ): 


i max min 
x.y v<zi 

subject to 

or equivalently, 

max 
u,v,t, 7 

subject to 


p(v) H G(X,V)p{ v ) 
tr ((S/(V) + ^Intn) *) 
tr (X) = E, X y 0 , V h 0 


t 

p{v) H G(U , V)p(v) > t, for VrGl 
tr ((S/(V) + ^ I Nt n ) U) = 1 
tr (17) = E'y, 7 > 0 


(13) 


(14) 


UhOVhO 
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where U = 'yX. Let g = [go,gi,--- ,gL~i] T with g t = 
J2k=i G(U, V)i+k,k, l = o, 1, • • • , L—l. One can also show 
that the constraint p{v) H Gp{v) > t in (14) is equivalent to 


f{v) = go - t + 2Re 

The optimization problem (14) is still non-convex and it 
includes infinitely many quadratic constraints as v £ X. To 
deal with this problem, we resort to an equivalent semi-definite 
programming (SDP) representation for the nonnegativity con¬ 
straint of the trigonometric polynomial in (15) based on [28, 
Theorem 3.4], which is quoted below as a lemma. 

Lemma 1: The trigonometric polynomial f(u>) = ho + 
2 Re{ Eti 1 hie~* ul } is non-negative over [a— f}, a+/3] (with 
0 < /? < 7 r) iff there exists an Lx L Hermitian matrix Z i V 0 
and an (L — 1) x (L — 1) Hermitian matrix Z 2 V 0 such that 


V 9ie- jlv 


> 0 


(15) 


h = F-l 1 (diag (FxZrff) + d © diag (F 2 Z 2 F 2 h )) (16) 

where h = [h 0 ,hi,--- ,h L - i] T , d = [d 0 ,di,--- ,dQ-i] T 
with d q = cos{2ttq/Q — a) — cos (/3), F\ = 
,/l-i] and F 2 = [/o, fi, • • • , fL- 2 ] where 
fi = [1, e ~ j2nl / Q , ■■■ , e -j27ri(Q-i)/Q] T with q > 2L - l. 

Based on Lemma 1, cyclic optimization [26] can then be 
performed to tackle problem (14) iteratively. To be specific, we 
perform the optimization with respect to U for some fixed V, 
and then conduct it with respect to V for fixed U, repeatedly. 
To this end, let a = vc, /? = A v and h = g te 1 in Lemma 
1, where e\ is an L x 1 vector with the first component being 
one and the others zero. For fixed V , the optimization with 
respect to U for (14) can be represented by the following SDP: 


max 

u,z 1 .z 2 ,t 

subject to 

< 


t 

g-t ei = F? (diag (Fi Z t F{ r ) 

+ d © diag {F 2 Z 2 F 2 )) 
tr ((£/(V) + ^P-IntN^U^J = 1 

u h 0 , z 1 y 0 , z 2 h 0 


(17) 


Let U* denote the optimal solution of U to (17). Then, the 
optimal solution of X is equal to EU */tr (U*). The SDP 
problem can be solved efficiently using the interior point 
methods in polynomial time [29]. In the simulations, we use 
the MATLAB toolbox CVX [30] to solve problem (17). 

Since the denominator of the objective function in (13) can 
also be expressed as tr ((S/(X) + In r n) V) with S/(X) = 
SfcLi INR kB(9k)XB(6k) H . the optimization with respect 
to V for fixed X can be cast as a similar SDP as below. 

max t 

v,z lt z 2 ,t 

subject to g — te 1 = F-[ r (diag (Fi Z-\ F-[ r ) 

< + dQdmg(F 2 Z 2 F 2 H )) (18) 

tr ((S/(X) + In r n) V) = 1 
V h 0 , Z 1 h 0 , Z 2 A 0 

By starting from a random initial point and cyclically 
solving (17) and (18) until the SINR improvement is negli¬ 
gible, the objective function value is non-decreasing and the 


convergence of the algorithm can be guaranteed [25]. The 
cyclic optimization converges to a point which is not only 
the local optimum, but also the global optimum along the X 
dimension and the V dimension separately [18]. To obtain 
a more accurate result, one can perform this procedure with 
a large number of random initializations and then select the 
best (X, V). In section IV, numerical examples show that the 
proposed algorithm is insensitive to initial values. 


B. Synthesis of s and w from X and V 

Let (X*, V*) denote the solution of (13) using the cyclic 
optimization. If both X * and V* are rank-one, the transmit 
waveform s * and receive filter w* can be obtained by the 
eigen-decomposition of X* = s*{s*) H and V* = w*(w*) H . 
In this case, the rank-one relaxation in (14) is tight and the 
solution is optimal. Otherwise, a suboptimal procedure can be 
adopted following a recently proposed algorithm in [25]. The 
basic idea of the algorithm is based on the fact that X(s,w) 
is a scaled version of the numerator p{v) H G{X , V)p(v), or 
tr (XA(d 0 ) H VA(d 0 )) equivalently. Then s* and w* should 
be designed to let |(tn*) H A(0 o )s*| 2 well approximate the 
shape of tr (X* A(9o) H V* A(9 q)) , while imposing constraint 
on the denominator. Interested readers can refer to [25] for 
detailed motivation. To make the letter self-contained, we shall 
present the synthesis algorithm for our problem in the sequel. 

Consider the value of tr (X*A(9q) h V*A(9q)) evaluated 
on DOAs (t?i, $ 2 , • • • , 9m} “uniformly distributed" on tt: 

c m = tr (X*A(9 m ) H V*A(9 m )) , m = 1,2, • ■ ■ ,M (19) 

Let T m = A(9 m ) H V*A(9 m ), QmQm = T m and define 
M auxiliary unit-norm vectors qi,q 2 ,--- ,qM■ Then, the 
synthesis of s can be formulated as 

{ min ELi II Qm,s - Jcffq m || 2 

s, Ql , ■ ,QM 

subject to s H ^E/(V*) + tr ^ ^ In t n^) s < (* (20) 

H<Zm|| = 1 , 1 <m< M 

where (* = tr((S/(V*) + tr ^ ^ In t n)X*). This problem 
can be solved using cyclic minimization. For a fixed .s, the 
solution to (20) is given by q m = yQ ^|| ,m = l,--- ,M. 
For fixed q rn ,rri = l,--- ,M, problem (20) reduces to a 
quadratically constrained quadratic program (QCQP) that can 
be solved by the CVX package [30]. The initial value of 
s can be chosen as the eigenvector of X* corresponding 
to the largest eigenvalue. Let s* denote the optimal solu¬ 
tion to (20), the optimal transmit waveform s* is given by 
s* = \EEs*/||s*||, considering the energy constraint on s. 

_ Analogously, let T m = A(9 m )X*A(9m) H , Q m Q% = 
T m , and 77* = tr ((S/(X*) + In r n) V*), the synthesis of 
w is similar to problem (20): 

! mm_ Eli \\Q m W - y/cflq 

m || 

w, Q lc ,Qm 

subject to w H (S/(X*) + In r n) w < if (21) 
llowll = 1, 1 <m< M 

which can also be solved using the cyclic minimization. 
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We also note that the randomized method [27] can also be 
used to obtain an approximate s* and w* in the non-rank-one 
case. Similar applications can be found in [20], [31], [32]. The 
synthesis algorithm based on the randomized method for our 
problem is shown in Algorithm 1. 


Algorithm 1 Synthesis algorithm based on randomized 
method_ 

Input: X * and V * 

Output: A randomized approximate solution s* and w * 

1: if rank(V*) = 1 then 

2: find w * via eigen-decomposition V * = 

3: else 

4: draw R random vectors Wj from the complex Gaussian 

distribution CAf(0, V*),j = 1,2 ,■■■ ,R 
5: calculate 

. w?A{9 0 )X*A(9 0 ) h Wj 

J e 0 eo Wj 1 Y,i(X*)wj + w^Wj ’ ^ ’ ’ 

where S/(X*) = £f =i INR, k B(0 k )X* B{0 k ) H . 

6: let to* = Wj mBX where 


Jmax — arg 


max A. 
i<j<H 


7: end if 

B: if rank(X*) = 1 then 

9: find s* via eigen-decomposition X* = s*(s*) H 

10 : else 

11: draw R random vectors s,; from the complex Gaussian 

distribution CJ\f(0 , X*), i = 1,2,--- , R 
12: calculate Si = and 


Q = min 


\(w*) H A(9 0 )si 


8 0 en (w*) H '£ii(si)w* + (in*) 


rk\H n 


* = !,•■• ,R 


where X)j(sj) = INR fc _B(0fc)sisf. B{B k ) H .. 

13: let .s' = s imax where 


*max = arg max Q. 

1 <i<R 

14: end if 


of {-2, -1,0,1,2} x {-60°, -50°, -40°, 40°, 60°, 70°}. The 
INR of all interferences is 30 dB. 

In Fig. 1, the output SINR as a function of 9q for the non- 
robust design and the proposed robust design are compared 
under four different parameters. For the non-robust design, 
the assumed a-prior target angle is set to be 9c and the 
optimization algorithm is based on the method presented in 
[19]. It is shown that the robust design improves the worst- 
case SINR performance significantly at the cost of peak-SINR 
degradation. For fixed A 9 and N, the superiority of robust 
design increases with the number of transmitters or receivers. 
In Fig. 2, we depict the beampattern P(9) = 


for parameter settings in Fig. 1(c) and Fig. 1(d) as an example. 
One can observe that both robust and non-robust design can 
produce nulls near the DOAs of interferences. From Fig. 1 and 
Fig. 2, we see that when A 9 is large enough relative to the 
beamwidth, the robust design can form a wide and flat beam 
over the uncertain space area to bring robustness. Both X* 
and V* are rank-one in this example. 



(a) 



■ robust 

- non-robust 

0 I—,—.—.—.—*—:—;—;—;— 
_10-8 -6 -4 -2 0 2 4 6 8 10 
0 O (degree) 

(b) 
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\ 
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« • 
i* 
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t 

i 

i 


-40. 

10 12 14 16 18 20 22 24 26 28 30 


e 0 (degree) 

(d) 


Fig. 1. Comparisons of the output SINR. A 9 = 10°. E = N = 20. 
SNR, = —15dB. (a) Nr = Nr = 4, 0c — 0°; (b) Nr = 4, Nr = 8, 
9c = 0°; (c) Nr = Nr = 8, 9q = 0°; (d) Nr = 8, Nr = 16, 9q = 20°. 


Prior results on the tightness of SDR [27], [33] show that for 
a separable SDP [27, eq. (28)] with P semi-definite variables 
and J constraints, there exists a rank-one optimal solution 
if J < P + 2. But this can not guarantee the existence of 
rank-one solution for our problem, since P = 3 and J = 
2L — 1 for problem (17) and (18) in the form of [27, eq. 
(28)]. Nevertheless, we emphasize that as with in [25], one 
can empirically observe that both X* and V* are rank-one 
for most of the random initializations as along as fl fl c = 0, 
where fl c denotes the set of all interferences angles. 

IV. Numerical Examples 

In this section, numerical examples are conducted to exam¬ 
ine the performance of the proposed method. In all examples, 
we assume that 30 interferences are present with the range and 
angle pair (i'k,0k) generated from all possible combinations 




Fig. 2. Comparisons of the beampattern P(9). (left) Nr = Nr = 8, 
9 C = 0°; (right) Nr = 8, N T = 16, 9 C = 20°. 

In Fig. 3, we plot the worst-case SINR versus the target 
angle uncertainty A 9. As expected, a wider range of target 
angle uncertainty leads to a worse SINR. The impact of A 9 on 
the worst-case SINR performance of non-robust design is more 
prominent, which suffers a sharp decline as A 9 increases. This 
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is due to the effect of the first null near the main lobe. In this 
example, both X* and V* are rank-one. 



A0 (degree) A9 (degree) 

Fig. 3. The worst-case output SINR versus the angle uncertainty. 6q = 0°. 
E = N = 20. SNR = -15dB. (left) Nr = 4 , N T = 8; (right) Nr = 
Nt = 8 . 




(a) (b) 




(c) (d) 


In Fig. 4, we investigate the effect of initial values on the 
cyclic optimization of X and V. We plot the worst-case SINR 
for the relaxed problem (14) under T = 50 different random 
initializations. Four different parameter settings are consid¬ 
ered. The cyclic optimization is stopped if either the increment 
of the worst-case SINR between two iterations is less than 
5 x 10 4 or the maximum number of iterations reaches. The 
maximum number of iterations of the cyclic optimization is set 
to 150. We can see that the worst-case SINRs under different 
initializations are very close. Let T = {t^\ t.^ 2 ' 1 , ■ ■ ■ , 
denote the worst-case SINRs from F random initializations. 
We define the following metric 


C = 


max(T) — rnin(T) 
mean(T) 


( 22 ) 


to evaluate the variation of T, where max(T), min(T) and 
mean(T) denote the maximum, minimum and mean value of 
T, respectively. The values of C for the four cases are equal 
to 0.016, 0.0177, 0.0263 and 0.016, respectively. One can see 
that in our problem, the cyclic optimization is quite insensitive 
to the initialization. 

In Fig. 5, we illustrate the performance of the synthesis 
algorithm in the non-rank-one case, which seldom happens in 
our experiments. In this example, the parameter settings are 
the same as in Fig. 4(c). Under a certain random initialization, 
the cyclic optimization provides a solution with rank(Tf') = 2 
and rank(V*) = 1. The receive filter w* is obtained based 
on eigen-decomposition, and the transmit waveform s * is 
obtained via the synthesis algorithm. The performance of 
synthesis algorithm based on solving problem (20) (denoted 
Method 1) and the algorithm based on randomized method 
(denoted Method 2) are compared. We plot their corresponding 
SINRs as a function of 0 o according to (2). For the Method 
1, the number of DOA samples M is set to 41 and the 
number of iterations to solve (20) is 50. For the Method 2, the 
number of random samples is set to be 1000. We also plot the 

SINR relax (0 o ) = tr((s V j(x , )+JWflJv)v ,) as a benchmark 
for comparison. We can observe that their SINR performance 
are very close, and both synthesis algorithms yield a good 
solution in the non-rank-one case. We can also see that the 


Fig. 4. The effect of initial values on cyclic optimization. E = NtN, 
SNR = 1/L, 6 C = 0° and A0 = 10°. (a) N R = N T = 4, N = 10; (b) 
N r = Nt = 4, N = 20; (c) N R = 4, N R = 8, N = 10. (d) N R = 4, 
N t = 8, N = 20. 


SINR curve of Method 1 matches well with SINR le i ax (0o)- The 
worst-case SINRs for Method 1, Method 2 and SINR re i ax (0o) 
are 18.478 dB, 18.394 dB 18.526 dB, respectively. 



e 0 (degree) 

Fig. 5. The SINR performance of the synthesis algorithm. 


V. Conclusions 

A method for angular-robust joint design of transmit wave¬ 
form and receive filter is proposed to maximize the worst- 
case SINR performance. The proposed method exhibits a con¬ 
siderable performance increment over the non-robust design 
via numerical examples. Future work will concentrate on the 
robust design with respect to the interferences uncertainty. 
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